Laser-assisted-autoionization dynamics of helium resonances 
with single attosecond pulses 



Wei-Chun Chu, 1 Song-Feng Zhao, 2 and C. D. Lin 1 

1 J. R. Macdonald Laboratory, Department of Physics, 
Kansas State University, Manhattan, Kansas 66506, USA 
2 Key Laboratory of Atomic and Molecular Physics and Functional Materials of Gansu Province, 
College of Physics and Electronic Engineering, Northwest Normal University, Lanzhou 730070, China 

(Dated: March 4, 2013) 

The strong coupling between two autoionizing states in helium is studied theoretically with the 
pump-probe scheme. An isolated 100-attosecond XUV pulse is used to excite helium near the 
resonance state in the presence of an intense infrared (IR) laser. The laser field introduces 
strong coupling between 2s2p( 1 P) and 2p 2 ( 1 5') states. The IR also can ionize helium from both 
autoionizing states. By changing the time delay between the XUV and the IR pulses, we investigated 
the photoelectron spectra near the two resonances. The results are used to explain the recent 
experiment by Gilbertson et al [Phys. Rev. Lett. 105, 263003 (2010)]. Using the same isolated 
attosecond pulse and a 540 nm laser, we also investigate the strong coupling between 2s2p( 1 P) 
and 2s 2 CS) by examining how the photoelectron spectra are modified vs the time delay and the 
possibility of observing Autler-Townes doublet in such experiments. 

PACS numbers: 32.80.Zb,32.80.Fb 



I. INTRODUCTION 

Quantum coherence is essential to the understanding 
and control of the dynamics in a quantum system. With 
the advances of experimental techniques and technolo- 
gies, extreme ultraviolet (XUV) attosecond pulse trains 
(APT) and single attosecond pulses (SAP) have been 
produced through the high-order harmonic generation 
(HHG) process, by exposing atoms to intense infrared 
lasers. Since the natural time scale of the electronic mo- 
tion of valence electrons in atoms and molecules is in 
the order of hundreds of attoseconds, these pulses can 
be used to manipulate or control electron dynamics at 
the attosecond level pj. For this purpose, pump-probe 
scheme is used experimentally where the dynamics of the 
system is investigated by varying the time delay between 
the pulses. However, existing XUV SAP and APT are 
too feeble to be used for XUV-pump-XUV-probe mea- 
surements. Instead, most of the experiments have been 
carried out using the weak XUV-pump with an intense 
IR-probe pulse. Rigorously speaking, the IR can be con- 
sidered as a probe only if it is applied after the XUV 
pump ends. In such case, the electron wave packet gen- 
erated by the APT or the SAP is probed. In most ex- 
periments, however, measurements are also made where 
the XUV and IR overlap in time. Since the IR is much 
stronger than the XUV, such measurements are better 
understood as laser-assisted photoionization. 

In XUV+IR experiments, the interaction of the 
XUV+IR field with the target atom is a nonlinear 
process. To treat such problems theoretically, the 
brute-force numerical solution of the time-dependent 
Schrodinger equation (TDSE) has been carried out by 
many groups where the target atom is treated within 
the single electron approximation. For APT+IR, Flo- 
quet approach 0] has also been employed as well. With 



much efforts, TDSE calculations have been performed 
for two-electron helium atom and two-electron hydrogen 
molecules and compared with measurements [3 5] . Still, 
simplified theoretical models are highly desirable in or- 
der to acquire a better understanding of the dynamics. 
These models, clearly, have to depend on the physical 
systems on hand. To be specific, we will focus on helium 
atom where many experiments have been carried out. 

The first type of XUV+IR experiments used APT 
where the photon energy runs from below to above the 
single ionization threshold of helium [1-Q3. Spectra of 
electrons or ions have been measured and calculations 
based on TDSE and Floquet theory reported 0, Q . In 
the limit of non-overlapping APT and IR, a simplified 
theory has been reported where the aim was to extract 
the electron wave packet generated by the APT [ll[ . Al- 
ternatively, photon spectra can also be measured fl2| . 
In this case, effect due to propagation in the medium 
has to be considered [l|. As the energy of the XUV 
photon increases beyond the singlc-ionization threshold, 
it enters a structureless spectral region (about 25 to 55 
cV). In this case, the addition of an intense IR to the 
XUV is to shift the momentum of the continuum elec- 
tron along the polarization direction of the IR laser, re- 
sulting in the "streaking" of photoelectrons where the 
energies depend on the vector potential of the IR laser at 
the time of the electron emission. A "streaking" theory 
based on the strong field approximation (SFA) [3] has 
been widely used. This theory also forms the basis of 
the frequency-resolved optical gating for complete recon- 
struction of attosecond bursts method [l5| for extracting 
the pulse duration of an SAP. At a still higher photon 
energy of 60 eV, the XUV alone will reach the spectral 
region where doubly-excited states of helium are located. 
In this paper, we will focus on this energy region and 
consider the dynamics of the doubly excited state gencr- 
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ated by the XUV pulse in the presence of an IR pulse at 
different delay times. 

Doubly excited states have been widelystudied using 
synchrotron radiations since the 1960's [HI [TtJ- These 
are autoionizing states where the lifetime of a state is 
typically of a few to tens of femtoseconds. These lifetimes 
are deduced from the measured spectral width; thus high- 
precision spectroscopic measurement is needed. With the 
emergence of SAP, the lifetime of an Auger resonance 
was first measured in the time domain and analyzed by 
the XUV+IR streaking theory in 2002 However, 
the evolution of the spectral shape of a resonance state 
has not been determined so far in the time-domain mea- 
surements, which in principle can be investigated with 
attosecond pulses. In fact, the autoionization theory was 
formulated in the energy domain byFano (l9| . This 
issue was addressed recently by us [20| and an experi- 
mental scheme was proposed. The scheme, however, re- 
quires an attosecond XUV-pump to create the resonance 
and an attosecond XUV-probe to project out the time- 
evolution of the resonance. Such measurement is not 
possible yet. Thus for the time domain measurements so 
far, doubly-excited states still have to be "probed" using 
an IR pulse. For an isolated Fano resonance generated 
in a combined XUV+IR pulse, the problem has been ex- 
amined by Wickenhauser et al (2~lj . A simpler model was 
proposed by Zhao et al [22j based on SFA. 

An XUV+IR experiment on helium doubly-excited 
states has been reported by Loh et al [23|, where a 30 
fs XUV pulse was used to excite the 2s2p( 1 P) doubly 
excited state in the presence of a 42 fs IR laser. The 
absorption spectra near the 2s2p resonance were mea- 
sured vs the time delay between the pulses. The result 
was then interpreted by a theory that includes the strong 
coupling between the 2s2p and 2p 2 { 1 S) states by the IR. 
In this experiment, the pulse duration of the XUV was 
longer than the lifetime of the 2s2p autoionizing state, 
i.e., the bandwidth of the light was narrower than the 
resonance width. Thus, the continuum electrons in a 
measurement did not have a meaningful "distribution" 
for such a narrow bandwidth. A single pump could not 
reveal any significant spectroscopic features; instead, se- 
quential measurements were made where the detuning of 
the XUV was controlled and scanned through an energy 
range. The scanned spectra were able to exhibit the au- 
toionization modified by the time-delayed IR, thus pro- 
viding the temporal information of the system partially. 
More recently, an 100 as SAP was used by Gilbertson et 
al [24[ to excite helium in the presence of a 9 fs IR pulse. 
The central energy of the excitation is at the 2s2p reso- 
nance, and the photoelectron spectra near this resonance 
vs time delay were reported. The result was interpreted 
based on the theory of Zhao et al [22j from which the 
decay lifetime of 2s2p was retrieved. However, the spec- 
tral shape of the resonance was not analyzed due to the 
limited resolution of 0.7 eV of the electron spectrometer. 
In contrast to Ref. (23|, the interpretation was made by 
the streaking model of Zhao et al [22[ only. Since the IR 



wavelength used in the two experiments are about the 
same, the role of strong coupling between the 2s2p and 
2p 2 states in this experiment should be addressed. This 
is the goal of the present paper. 

The present analysis limits the Hilbcrt space to includ- 
ing only the two autoionizing states, 2s2p and 2p 2 , and 
the ground state. While the XUV can populate other 
higher doubly excited states, they are weaker in the spec- 
trum and are easily ionized by the IR. Governed by the 
detuning of the 9 fs IR, only the 2s2p and 2p 2 states 
are strongly coupled. The dynamics of such a system 
is then studied with a theory generalized from the stan- 
dard three-level system formulated earlier [25l - [28j . To 
be specific, we develop a model for the IR-dressed au- 
toionization of the 2s2p resonance in helium excited by 
an SAP, where the IR strongly couples the 2s2p and 2p 2 
resonances and also ionizes both doubly-excited states. 
Both effects by the IR have to be included in order to ex- 
plain the observed electron spectra reported in Ref. 24 1 . 
After achieving good agreement with the experiment, we 
further investigate the case where 2s2p and 2s 2 ( x 5) are 
coupled by changing the laser wavelength to 540 nm. The 
latter state has a higher binding energy such that it is 
not ionized by the laser. Such measurements would be 
closer to the standard electromagnetically induced trans- 
parency (EIT) experiment [1^, y(| . We will look for the 
presence of the Autler-Townes doublet [3l[ which are rou- 
tinely generated by three-level systems with long pulses. 

In Sec.HH the model system is defined and the method 
is introduced. In Sec. lIIIl we show and analyze the results 
of the two cases, with laser wavelengths A^ = 780 and 
540 nm, respectively For Al = 780 nm, our calculation 
is compared with the available experiment [24[ and with 
the streaking model based on the SFA model of Zhao et 
al [22j |. For Xl — 540 nm, the results are analyzed, and 
the possible experimental realization are discussed. In 
Sec. IIV1 we give the conclusion. Atomic units (a.u.) are 
used in Sec. [TTJ In the rest of the paper, electron Volts 
(eV) and femtoseconds (fs) are used for energy and time, 
respectively, unless otherwise specified. 



II. THEORY 
A. General description of the model system 

Consider an atomic system with ground state \g) and 
two doubly excited states \a) and \b). The dipole transi- 
tion is allowed between \g) and \a) and between \a) and 
| b). Suppose their energy levels and the external fields 
are arranged in a way such that the XUV is near reso- 
nance between \g) and \a) and the laser is near resonance 
between \a) and \b). The total time-dependent wave func- 
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tion of the system can be approximately written as 
Mt)) = e- iE ° t c g (t)\g) 



-iE x t 



-iE L t 



d a (t)\a)+ I d El {t)\E 1 )dE 1 
(k(t)\b) + ( d E2 (t)\E2)dE 2 



(1) 



where E g is the ground state energy, Ex = E g 
and El = E g + U)x +&l are the "pumped energies" with 
respect to the photon energies cox and ul, where X is for 
XUV and L is for laser; \E\) and \E 2 ) are the continua 
with respect to \a) and to \b). The fast oscillating parts 
in the wave function have been factored out, and the c(t) 
and d{t) functions are assumed slowly- varying with time. 
By convention, we use symbol c(t) for coefficients if the 
expansion is with respect to eigenstates, and symbol d(t) 
is used when the expansion is in terms of configurations 
(not eigenstates). The total Hamiltonian of the atomic 
system in the external holds is 



H(t)=H A + H x (t)+H L (t), 



(2) 



where Ha is the atomic Hamiltonian, and H x (t) and 
-ffi(i) represent the interactions of the XUV and laser 
fields with the atom, respectively. In the photon energy 
range in consideration, the interaction with the field is 
given by the electric dipole transition. Equation [2] pro- 
vides the total Hamiltonian to solve the time-dependent 
Schrodinger equation for the wave function in Eq. [TJ i.e., 



mm))- 
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The fields are assumed linearly polarized in the same di- 
rection. The electric fields arc in the form of 



E{t) = 2F(t) cos(o;i) 

= F{t) (e iujt + e- lut ) 



(4) 



where 2F{t) is the pulse envelope. The present model 
works for arbitrary pulse envelopes; however, we limit 
the envelope to cosine-square shape in the calculation, 
i.e., 



F{t) = F cos 2 



t-tn 



7TT 7TT . . 

for - — < t - t Q < — (5) 



and F(t) = anywhere else, where to is the peak time. 
For such a pulse, the pulse duration is 7rr/2.75, and Fo 
is related to the peak intensity Io by Io = 4Fq . In this 
paper, the pulse duration is defined by the full width at 
half maximum (FWHM) of the intensity envelope. The 
time delay between the two pulses is define as the time 
of laser peak subtracted by the time of XUV peak, so it 
is positive if the peak of the XUV appears before the IR. 
By convention, in this paper, the XUV peak and the IR 
peak are placed at t — and t = to respectively, so the 
time delay is to- 




FIG. 1: (Color online) Diagram for off-diagonal terms of the 
Hamiltonian in the configuration basis. The double-arrows 
represent the transitions between the basis vectors; the red 
ones represent the dipole transitions, the green ones represent 
the CIs responsible for autoionization, and the gray dashed 
ones are neglected in the model. 



Note that the basis functions for the autoionizing 
states in Eq. Q] are not the eigenstates of the atomic 
Hamiltonian Ha ■ In the basis space of Eq. [l] the diag- 
onal terms of the Hamiltonian are just the energies E g , 
E ai E\, Eb, and E 2 . The off-diagonal terms are 



(Ei\H(t)\a 
{E 2 \H(t)\b 

{a\H(t)\g 
(Ei\H(t)\g 

(b\H{t)\a 
(E 2 \H(t)\a 



V a 
V b , 



-iujxt 



-D ag F x {t)e 
-Di g F x (t)e- iuxt 
-DbaFL^e-' 1 ^ 
-D 2a F L {t)e- iULt 



(6) 
(7) 
(8) 
(9) 
(10) 
(11) 



where the D are the dipole matrix elements and V the 
transition amplitudes of configuration interaction (CI). 
These terms arc schematically plotted in Fig. [TJ 

For the process in concern, additional simplifications 
are made. Helium atom is taken as a prototype, where 
Iff) is lsY-S"), I a) and \b) are 2s2p( 1 P) and 2p 2 ( 1 S), 
and \Ei) and \E 2 ) are lsep( 1 P) and lsesQ-S), respec- 
tively. Referring to Fig. [TJ since the IR coupling is a one- 
electron dipole operator, the coupling between 2p 2 ( 1 S') 
and lsep( 1 P) is zero to the first order. The coupling 
between the two continuum states by the dipole opera- 
tor can also be neglected since the IR is not absorbed 
by the continuum electron. We also use the rotating 
wave approximation such that only the resonant transi- 
tions are considered. Furthermore, the matrix elements 
involving continuum states \Ei) and \E 2 ) are assumed 
energy- independent, which means D and the Fano pa- 
rameters, r and q, are constant values estimated at the 
resonances. This is a good approximation when the reso- 
nance energy is high enough above the threshold relative 
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to its width so that the continuum only varies slightly 
across the resonance. In Eqs. IMTTl the basis functions 
are real by convention (the continua are standing waves) 
so that all the D and V are real. In this model, all the 
atomic parameters (D, E, T, and q) are taken from ex- 
periments or from calculations in the literature. These 
parameters are determined by the atomic structure which 
is irrelevant to the development of our model. 



B. Three-state model 



that the model described above retrieves only the bound- 
state part of the total wave function. It has been ap- 
plied successfully for long pulses where the photoelec- 
tron distribution in each single measurement is disre- 
garded [H, HH-HH. In the present study, however, a 
broadband XUV pumps the electrons to the continuum 
and to the bound state alike; missing the continuum co- 
efficients means missing the essential information of the 
total wave function. In the following, wc aim to extend 
the model to recover the dynamics of the complete sys- 
tem for short-pulse cases. 



Under the approximations outlined above, a system 
of coupled equations for all the coefficients appearing in 
Eq. [1] arc obtained. Considering the conservation of the 
total probability of the wave packet, the continuum-state 
coefficients change with time much more slowly than the 
bound-state coefficients. Thus, d El (t) and d E2 (t) are adi- 
abatically eliminated by assuming their time-derivatives 
are zeros [Ul^. This allows us to reduce the calculation 
to include only the bound states. Solving the coupled 
equations is then numerically feasible. Now we have 



•7 9 (t) 



2 " 9 

id a (t) = \ a F x {t)c g {t) - 
+ X b F L (t)d b (t) 



c g {t) + X a Fx(t)d a (t) 



(12) 



5 X + i 



dait) 



(13) 



id b [t) = X b F L {t)d b (t) - \S X +S L + <k(t), (14) 

where A a = -D ag (l - i/q a ) and A 5 = -D ba (l - i/q b ) 
are the newly defined complex dipole matrix elements 
which include the route through the continua. The laser- 
induced broadening r y g (t) and 7 a (£) are defined by 



lg {t) = 2-K\D lg F x {t)Y 
7o (t) = 2n\D 2a F L (t)\ 2 



(15) 
(16) 

E q — E a and 



the detuning of the pulses are 5x = 
5l = + E a — E b , and the Fano parameters are r a 
2n\V a \ 2 , T b = 2v\V b \ 2 , q a = D ag /{TrV a D lg ), and q b 
D ba /{^V b D 2a ). The AC Stark shifts vanish, i.e., 



\DigFx(t)f 
Ex 



E 



dEx = P 



x 



\D2aF L (t)\ 

El — E-2 



■dE 2 = (17) 



because D\ g and D 2a are constant of energy, where P 
means the principal part. 

With any given set of atomic parameters, field pa- 
rameters, and initial conditions, Eqs. I12M14I uniquely de- 
termine the bound-state part of the total wave func- 
tion regardless of the continua. In numerical calcula- 
tions, we propagate these coefficients using the Runge- 
Kutta method, with the initial conditions c g (ti) = 1 and 
d a (ti) = d b (ti) = 0, i.e., the system is in the ground 
state at initial time U before the pulses arrive. Note 



C. Continuum states 

The original coupled equations for the continuum 
states before the adiabatic elimination are 

id El {t) = -D lg F x (t)c g (t) + V a d a (t) + {Ex - E x )d El (t) 

(18) 

id E M = -D 2a F L (t)d a (t) + V b d b (t) + (E 2 - E L )d E2 {t). 

(19) 

Assuming d El (t) = d E2 (t) = in Eqs. [T8]fT9l the approx- 
imate d El (t) and d E2 (t) have singularities at Ex = E x 
and at E 2 = El, respectively. When plugged into the 
bound-state coupled equations, these singularities are 
handled properly with contour integrations. After solv- 
ing Eq. [T21IT41 the bound-state coefficients c g (t), d a (i), 
and d b (t) are known functions of time. We then return 
to the full version of the continua in Eq. [T8llT9l to retrieve 
d El (t) and d E2 (t), again by the Runge-Kutta propagation 
over time, with the initial conditions d El (ti) = d Ez (ti) — 
0. The new d El (t) and d E2 (t) functions generated by 
Eq. [T8HT91 are the next iteration and better solutions to 
the previous ones. 

After the pulses are over, the system may still change 
due to autoionization. For an arbitrary set of bound- and 
continuum-state coefficients, the photoelectrons evolve to 
a final energy distribution after about 10 times the de- 
cay lifetime or so [20| . To simulate this measurable pho- 
toelectron spectrum, we have to project the total wave 
packet that has propagated for a long time onto the par- 
tial waves of interest; equivalently, taking the |a)-|£i) 
resonance for example, the spectrum is proportional to 
Mbi(^/)| 2 as if ~ * 00 ■ The probability density of the 
continuum wave function, P{Ex) = \d El {tf)\ 2 , is called 
"photoelectron profile" or simply "profile" hereafter in 
this paper. Numerically, tt is taken when Id^^i)! 2 does 
not change anymore. In the case concerning the 2s2p( 1 P) 
resonance in helium, the decay lifetime is 17 fs, so the fi- 
nal spectrum can be safely taken at about tf = 200 fs. 
Note that P{Ex) has the dimension of probability per 
unit of energy, and the integral of P{Ex) over energy 
represents the total probability of the system being in 
the I -Ei) continuum. The conservation of total proba- 
bility can be used to numerically check the convergence 
of calculation. The continuum states calculated this way 
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have some advantages. First, each of them can be turned 
on separately once the bound-state calculation is done. 
Second, the energy range and energy mesh can be ar- 
bitrarily chosen without lowering the accuracy because 
d.E 1 {t) and dE 2 (t) are calculated independently for each 
energy point. 

In the present work, we only concern ourselves with 
the final photoclcctrons reaching the detector. However, 
if one wants to know the evolution of electrons in real 
time before the end of decay, the retrieval of ^^(i)! 2 
and \dE 2 (t)\ 2 as functions of energy and time will be the 
answer, instead of P(E\) and P{E2) as just functions 
of energy. The measurement of this short-time behavior 
was proposed using an additional high-energy short pulse 
to ionize the inner electron 1201. 



D. Eigenstates 

The atomic eigenstates near a resonance are solved in 
terms of the corresponding bound state and background 
continuum by Fano's theory (l9| . This leads us to express 
the total wave function in eigenstate basis in the general 
form of 



-iExt 



cf{t)\^)dE 



+ e- iE -* I c$m%))dE>, 



(20) 



where the superscripts (a) and (b) indicate the |a}-|-Ei) 
pair and the \b)-\E2) pair, respectively. The eigenstates 
we refer to are the eigenstates of atomic Hamiltonian Ha, 
so they are stationary only in the absence of field. In the 
following, we will discuss only the subspace spanned by 
| a) and \E\) in details, while the same principles apply 
to 1 6) and -E^)- In the meantime the superscript (a) is 
omitted. 

The solutions for an eigenstate of energy E in config- 
uration basis is 

He) = His\a) + J vee x \Ei)dE\, (21) 



where the coefficients are 
sin 9e 



[IE = 

vee x = 



sin 9e 
tt(E - Ei) 



(22) 

-cos6 E 6(E - E±), (23) 



where 



9e = — tan 



i ( r a /2 



E-E a 



(24) 



The conversion between Eq.Q]and Eq. [20] means that for 
the | a) resonance, 



ce 



(t)\ifi E )dE = d a (t)\a) + J d El {t)\E 1 )dE 1 . (25) 



Starting with Eq. [21] and Eq. 

ce(t) coefficients are 



with some algebra, the 



c E {t) 



sin U e 



d a (t) - (cosflg - ism0 E )d El (t)\ _ 



E 1= E' 



(26) 

which enables us to calculate the detected photoclcctron 
spectrum P(E) = ce(i/)| 2 as tf — » oo. This spectrum 
defined here will reach the same result as P{Ei) defined 
in Sec. Ill CI but just from a different calculation proce- 
dure. 

The profile |ce(£)| 2 evolves differently from | (*) | 2 ■ 
As pointed out in Sec. Ill CI Ide^i)! 2 evolves until the 
end of the decay, and |cE(t)| 2 evolves until the field is 
over. In a typical case where the field ends much earlier 
than the end of the decay, it is more efficient to calculate 
the final |ce(£/)| 2 than the final \d El (i/)| 2 - Moreover, 
in Eq. 1261 for a point in time and in energy (t,E), the 
coefficient ce(t) is calculated by a simple algebra with 
the d a {t) and dE 1 {t)\E 1 =E of just the same (t,E); there 
is no integral over energy nor propagation over time to 
carry out cg:(i) once d a (t) and dE 1 (t)\E 1 =E are known. 
In other words, including ce(£) in the calculation requires 
very little extra efforts. Summing up these facts, for the 
purpose of retrieving the electron spectra, it is advanta- 
geous to adopt P(E) = |ce(i/)| 2 instead of keeping the 
formofP(£i) = \dE 1 (tf)\ 2 - We have come to the conclu- 
sion that for the present system and setup, the optimal 
workflow is in the consecutive order of bound states, con- 
tinuum states, and eigenstates, all of which calculated in 
the physical time until the field vanishes. 

For the subspace spanned by \b) and \E<i), the eigen- 
state coefficients are constructed by db(t) and d_E 2 (t) in 
the same way but with different parameters. Now both 



(i) and 4 6) (t) are obtained, and we can fully build 
the total wave function in eigenstates basis as shown in 
Eq.HDl 



III. RESULTS AND DISCUSSIONS 
A. Laser wavelength of 780 nm 

1. Calculation and analysis 

We first apply the model to the case where the laser 
couples 2s2p( 1 P) and 2p 2 ( 1 S') in helium with the exper- 
imental setup reported by Gilbcrtson et al |24j. The 
two doubly excited states autoionize to lsep( 1 P) and 
Isesi^S) respectively. The measured quantity is the elec- 
tron spectrum, corresponding to P(E) = |cE(t/)| 2 in the 
model introduced in Sec. Ill Dl where tf is taken after the 
field vanishes. For the XUV pulse, the photon energy is 
60 eV and the duration is 100 as; the bandwidth is 20 eV, 
which is much wider than the width of 2s2p, and viewed 
as a flat background in spectrum near the resonance. The 
XUV transition from the ls 2 ( 1 S') ground state to 2s2p is 
nearly resonant, and we conveniently set Sx = 0. In 
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principle, such pump pulse can initiate other 2snp, 2pns 
and 2pnd states which should all be included in the total 
wave function. (For more accurate description of dou- 
bly excited states, see Lin [32j]) Nonetheless, those higher 
states have longer lifetimes, exhibiting less dynamics in 
the time scale in our scheme; furthermore, their reso- 
nance widths are narrow and difficult to measure with 
typical electron spectrometers. Thus, we treat 2s2p as 
the only \a) state in Eq. [TJ while Isep is \Ei). The Fano 
parameters of the 2s2p resonance are experimental val- 
ues taken from the literature [HI, where E a = 60.15 eV, 
T a = 37 meV, and q a = -2.75. The XUV is weak and 
in the linear regime so that its intensity does not affect 
the spectra other than an overall factor. The peak XUV 
intensity is assumed Ix = 10 10 W/cm 2 . 

The laser pulse is a 9 fs IR pulse with wavelength Xl = 
780 nm and peak intensity I L = 7 x 10 11 W/cm 2 . With 
1.6 eV photon energy, the IR field couples 2s2p and 2p 2 
at 62.06 eV, where Sl = —0.4 eV. The atomic parameters 
T h = 5.9 meV 0] and D ba = 2.17 a.u. [H are taken from 
literature. The bound-free dipole matrix element T>i a 
responsible for the transition from 2s2p to lses is very 
small because it is a second order (satellite) transition 
and requires electron correlation, while D ba between 2s2p 
and 2p 2 is a first-order transition via (2s\D\2p). The 
parameter q b , representing the ratio of D ba to Z?2a, is 
thus very large; it is set q b = 1000. Note that when 
\q\ 3> 1, the Fano shape appears as a symmetric peak, 
and the sign of q is insignificant. 

In order to understand the dynamics controlled by the 
laser, we evaluate the generalized Rabi frequency defined 
by 



(27) 



When the detuning is large, the Rabi frequency is higher; 
however, the amplitude of the oscillation is lower, i.e., 
the population does not fully swing to the other coupled 
state. In our system and setup, SI = 0.43 eV between 
2s2p and 2p 2 at the laser peak, corresponding to 9.6 fs 
period. For the two doubly excited states, the binding 
energies of 2s2p and 2p 2 are 5.3 and 3.3 eV, which are 
low enough such that their ionizations are quick, espe- 
cially the latter one. We calculate the ionization rates 
for both states using the model developed by Perelomov, 
Popov, and Terent'ev 35|, referred by the PPT model 
hereafter, with the high-intensity correction introduced 
in (36j . The empirical parameters in our PPT calcu- 
lation are obtained by fitting our result to the experi- 
ment, which will be discussed in Sec. MI A 21 The calcu- 
lated peak rates for 2s2p and 2p 2 , in units of energies, 
are 5.4 meV and 0.46 eV, comparable to their resonance 
widths 37 meV and 5.9 meV, respectively. To incorporate 
these ionization rates to our coupling model in Eqs. 1121 
ITU the widths are broadened by T' a b (t) = T a , b + W a , b (t) 
where W a (t) and W b (t) are the time-dependent ioniza- 
tion rates for 2s2p and 2p 2 respectively. Note that for 
the dynamics we have considered so far, on one hand, 
the coupling and the ionization of resonances exist only 
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FIG. 2: (Color online) Probabilities of the (a) 2s2p and (b) 
2p 2 bound states with time, for Xl = 780 nm and various to- 
The case without laser is plotted in the gray solid curve in 
(a). 



in the presence of laser field; on the other hand, the au- 
toionization processes are determined exclusively by the 
atomic structure, whose time scale can not be changed 
externally. 

The probabilities of the 2s2p and 2p 2 bound states 
are presented in Fig. [5] as they propagate in time. Al- 
though the bound states are not directly measured in 
this pump-probe scheme, the analysis therein helps re- 
veal the physics behind the whole time-dependent pro- 
cess. We will analyze the bound-state propagation and 
the resonance profiles in the following at the same time. 

The photoelectron profiles of the 2s2p resonance for 
to = —10 to 50 fs are shown in Fig. [3] If the laser field 
is absent, the dynamics of the system after the pump 
will be nothing more than the autoionization of 2s2p, 
and the detector will see an original Fano line-shape as 
seen in the absorption spectrum in synchrotron radia- 
tion experiment. When a laser is added very early, e.g., 
to = —10 fs, the spectrum only changes insubstantially 
because the laser is already gone by the time the XUV 
pumps the system. The laser has essentially no influ- 
ence on the dynamics; this is seem at the left end of the 
spectrogram. If the laser is shifted after the XUV, the 
laser will start pumping the system from 2s2p to 2p 2 . 
The period of Rabi oscillation is 9.6 fs, close to the laser 
duration 9 fs. The ionization of 2p 2 has the time scale 
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FIG. 3: (Color online) Photoelectron profile of the 2s2p reso- 
nance for Xl = 780 nm. (a) The spectrogram. The cyan bar 
indicates the range of delay where the two pulses overlap, (b) 
The spectra for various to, and the case without laser. 



= 1.4 fs estimated at the laser peak. Such a short 
time indicates that 2p 2 is very quickly depleted in the 
presence of the laser. For to = to 5 fs, the laser is 
mainly at the beginning of decay. Most population is 
brought from 2s2p to 2p 2 before 2s2p decays, and being 
ionized from 2p 2 without returning to 2s2p. Figure [2a) 
shows that for to = 2 fs, there is a quick drop at around 
t = 4 f s which takes away about 40% of the original pop- 
ulation if compared to the "no laser" level. The amount 
of the 2s2p doubly excited state is greatly reduced, gen- 
crating much less photoclcctrons and creating a signifi- 
cantly lower peak in the spectrum, as seen in Fig. EJa) . 
Figure [3jb) shows the change of the resonance profile vs 
the time delays in this overlapping region. If the laser 
moves further positively, the longer lags between the two 
pulses gives more time for 2s2p to decay before the cou- 
pling kicks in. For to 3> 17 fs, the laser is late enough so 
that 2s2p decays completely without being interrupted, 
i.e., the original Fano shape is fully restored, and the 
laser has nothing to pump anymore. The recovery of the 
resonance profile vs time delay is shown in Fig. (3jc). Ul- 
timately, because the laser only removes bound electrons, 
changing the delay time traces out the decay process. 

Figure @] shows the 2p 2 resonance profiles. The signal 
is large in the range from to = — 10 fs to and low for 
to > 15 fs. Starting from the left end of the spectrogram, 



FIG. 4: (Color online) As Fig. [3] but for the 2p 2 resonance. 



for to = —4 fs, the tail part of the laser pulse is involved 
in the dynamics; counting only this involved part of the 
laser, the strength is low and the duration is short, and 
the Rabi oscillation is slow. As a result, the maximum 
population pumped to 2p 2 is "trapped" without return- 
ing back to 2s2p. At the same time, the ionization rate 
by this laser "tail" is low and not influential. The high 
population of the 2p 2 state is shown in Fig. E£b). At 
large time, it decays into the biggest photoelectron pro- 
file shown in Fig. 2Jb). Moving on to to = 2 fs, the laser 
strikes 2s2p mainly at the beginning of decay. As shown 
in Fig.[2jb), the full strength of laser depletes the popula- 
tion almost completely; at the end of laser, the 2p 2 bound 
state is almost empty, and there are hardly any electrons 
to autoionize. Consequently, the 2p 2 profile in Fig. [1] is 
greatly depressed, forming a valley in the spectrogram. 

Note that the 2s2p and the 2p 2 resonances are in dif- 
ferent symmetries. Their momentum (or angular) dis- 
tributions are different, but their energy ranges over- 
lap because of the broad bandwidth of the SAP. If the 
momentum spectrum is measured, one can distinguish 
the contributions from the two resonances and sepa- 
rate out the spectra. However, if only the energy spec- 
trum is measured, the two contributions are not separa- 
ble near 37.5 eV, where the 2s2p resonance lies on the 
lsep background. The Isep background at 37.5 eV is 
7.5 x 10 -9 cV _1 , which is about the same as the peak 
value of the 2p 2 profile shown in Fig. 2] The 2p 2 res- 
onance profile in the total energy spectrum will appear 
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FIG. 5: (Color online) Photoelectron profiles for Xl = 780 nm 
and various II, of (a) the 2s2p resonance for to = 30 fs, and 
(b) the 2p 2 resonance for to = 5 fs. 



flatter than how it would look without the background. 

The mechanism of the dressing laser is further ana- 
lyzed when the calculations are done for various laser 
intensities. The curves in Fig. [5ja) and (b) represent the 
spectra for fixed delays but different intensities, where 
the medium one (II = 0.7 TW/cm 2 ) repeats the experi- 
mental condition. The delays 30 and 5 fs are chosen for 
2s2p and 2p 2 respectively for adequate signal strengths. 
For 2s2p, with the increase of Il 1 the profile gradually 
depresses and forms interference patterns. It suggests 
that if the laser at about 30 fs is strong enough, it will 
suddenly deplete the bound states and halt the decay pro- 
cess. The resultant profile is then constructed by the elec- 
trons autoionized before 30 fs. Fig. [5ja) also shows that 
the resonance peak moves gradually to the low-energy 
side with increasing IR coupling strength. For the 2p 2 
resonance shown in Fig. [5Jb) , increasing I L mainly pop- 
ulates more 2p 2 and makes the profile larger. 



2. Comparison with experiment 

The experiment by Gilbertson et al [24| reported the 
spectrogram in the energy range 33-46 eV, enclosing both 
resonances in our concern. However, the energy resolu- 
tion is insufficient for the very narrow 2p 2 resonance with 
Tb = 5.9 mcV, making its transient spectra invisible. The 
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FIG. 6: (Color online) The normalized peak value of the 2s2p 
profile as a function of delay for Xl ~ 780 nm. The value for 
to = — 10 fs is set as 1. The present calculations are shown in 
three different settings where the ionization rates W a {t) and 
Wb{t) are included or excluded. The experimental curve [24| 
is shifted in to by 4.5 fs to fit the full calculation shown by 
the red solid curve. The result by the SFA model [52] is 
normalized such that it fits the experimental data the most. 
The cyan bar indicates the delay range where the two pulses 
overlap. 



most dominant feature in the spectrogram is the varia- 
tion in peak height of 2s2p resonance, where the decay 
lifetime 17 fs can be extracted. In Fig. [51 the normalized 
peak value is shown as a function of to ■ The peak value 
for the most negative delay, which represents the "XUV 
only" case, is normalized to 1. The present calculation 
is displayed in three different settings, where the model 
excludes the laser ionization, includes the laser ioniza- 
tion from 2s2p, and includes the laser ionization from 
both states. While the experiment did not know the ab- 
solute times of the pulses thus unable to determine the 
absolute delay, the experimental data is shifted by 4.5 fs 
to fit our calculation. In Fig. [51 the calculation result 
without laser ionization exhibits the similar decay fea- 
ture found in the experiment, but the absolute value is 
too high. By adding the laser ionization, we are able to 
adjust the parameters such that the calculation agrees 
with the experiment quantitatively. The empirical ion- 
ization parameters are Ci = 0.5 and a = for 2s2p and 
Ci = 0.5 and a = 11 for 2p 2 , where Ci is defined in 
the PPT model (HJ and a is defined in the correction 
term [36j . 

The comparison shown above indicates that the pres- 
ence of the 2p 2 state is responsible for the significant de- 
pression of the Fano peak, where both the coupling and 
the IR ionization take credit. The agreement between 
our model and the experiment is very good when both 
effects are taken into account. However, there are irreg- 
ular oscillations in the experimental data; their period is 
typically between 5 and 10 fs, but mostly unpredictable. 
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Their appearance is unexplained so far, but could be due 
to experimental artifact. 

The experimental report included a simulation using 
the "streaking" model (within the SFA) for isolated au- 
toionizing states [22| , with the inclusion of ionization due 
to the laser calculated within the PPT model [35[ , as de- 
scribed by Eq. 1 in Ref. [H[ . In the SFA model, the 
ionization of 2s2p and the acceleration of the scattered 
electrons are treated separately, and the laser coupling 
to other bound or resonance states is totally disregarded. 
The common mechanism between the model therein and 
our model is that the laser ionization is treated by the 
PPT calculation. In order to have a clear view on the two 
approaches, we have reproduced the angular-differential 
electron spectrogram at the polarization direction, where 
the measurement was done, using the SFA model with the 
same parameters taken in our model. The result is shown 
in Fig. [7] The corresponding resonance peak is plotted in 
Fig. [5] which is normalized to fit the overall experimental 
data. The SFA results in Fig. [6] and in Fig. [7] are similar 
to our calculation for i > 10 fs; the resonance peak goes 
through the same depression before the gradual revival 
along to. However, for to < 10 fs, the strong streaking 
peaks carrying the laser period 2.6 fs are missing both 
in the experiment and in our model. This suggests that 
with the inclusion of coupling and ionization of 2p 2 , the 
present model takes charge of the primary effects seen in 
the time-delayed measurement. We comment that even 
though both the SFA model and the present one account 
for the ionization by the IR laser, it is the ionization of 
2p 2 that is mostly responsible for the great reduction in 
the resonance strength observed in the experiment. This 
state was not included in the SFA model. 

The experiment also measured the spectrograms with 
several other laser intensities, and reported the depth 
of the "dip" that is seen in Fig. [6js curve against 1^. 
The experimental and theoretical dip values are shown in 



Fig. IS] As seen in the figure, increasing 1^ will aggravate 
the depression of the resonance. The experimental dip is 
also deeper than the theoretical one for high II, i.e., the 
resonance peak is more depleted than what the present 
theory predicts. This is because as II increases, the non- 
linear interaction of the IR opens up more excitation and 
ionization pathways that deplete the 2s2p part (or the 
bound part) of the resonance. The depletion will lead 
to fewer autoionization events and thus weaker electron 
signal near the resonance. Note that the discrepancy be- 
tween the present model and experimental data in Fig. 8 
grows nonlinearly with II for II > 4.5 x 10 11 W/cm 2 . 



B. Laser wavelength of 540 nm 

Using long pulses and laser spectroscopy, Autler- 
Townes doublet [3(1 [HJ has been observed when two 
bound states are strongly coupled by a dressing laser 
field. For the two autoionizing states and the ultrashort 
laser in the present study, can we observe anything re- 
sembling the Autler-Townes doublet in the spectra? In 
Sec. IIII Al we have shown that for the 780 nm IR, since 
the laser detuning is large and the laser ionization is sig- 
nificant, the evidence of Autler-Townes doublets is not 
prevalent except that the small shift of the 2s2p peak at 
higher coupling intensity in Fig. [SJa) gives a feeble hint. 
Now, we intentionally tune the laser to Al = 540 nm 
(mean photon energy is 2.30 eV) so that it is in reso- 
nance between 2s2p and the lower-energy state, 2s 2 ( 1 S l ) 
at 57.85 cV. In this setup, first of all, the laser detuning 
is negligible, making the Rabi oscillation the strongest. 
Secondly, the binding energy of the \b) state changes from 
3.34 eV for 2p 2 to 7.55 eV for 2s 2 , which effectively shuts 
down its laser ionization. As a consequence, the Rabi 
flopping dominates, and other complications are mini- 
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mizcd. Below we examine whether we can observe the 
Autler-Townes doublet for such a "three-level system" 
where the two fast-decaying " levels" are strongly coupled 
by a short pulse. 

For the parameters in the model, the 2s 2 resonance 
width (r&) of 0.125 eV, or the lifetime of 5.3 fs, is taken 
from the earlier calculation [37]]. The dipole transition 
and the g-parameter are assumed to be the same as those 
for 2p 2 , i.e., D^a = 2.17 a.u. and qb = 1000, because the 
first-order transition in D^ a is again (2s\D\2p), and D 2a 
has no first-order term. The PPT rate for 2s 2 at the laser 
peak is in the order of 10 -5 eV which is negligible for our 
purpose. The Rabi frequency between 2s2p and 2s 2 is 
0.27 eV, corresponding to the period of 15 fs. 

The evolution of the bound states is shown in Fig. |9] 
An obvious difference from the Al = 780 nm case is that 
the laser coupling now is strong enough to completely 
deplete 2s2p so that its population almost touches zero 
before bouncing back, as shown in Figs. [5Jb-d). How- 
ever, when it revives, the amount brought back by the 
oscillation is only less than 10% of what has been re- 
moved. For the 2s 2 state for to = 5, 10, and 20 fs, the 
overall population decreases with the delay because the 
2s2p decays with time and the laser pumps less and less 
electrons from 2s2p to 2s 2 . 

Now we turn to the electron spectra of the 2s2p reso- 
nance in Fig. [TD] When the delay passes the origin and 
becomes positive, the resonance profile flips horizontally, 
looking like the mirror image of the "no laser" spectrum 
other than an overall reduction in height. This flipped 
peak is seen in the to — 5 fs curve in Fig. HUT c). To 
understand this pattern, in Fig. IHJb), the sharp drop of 



FIG. 10: (Color online) As Fig. H but for X L = 540 nm. The 
gray solid curve represents the "no laser" case. 



2s2p from to 7 fs suggests that most electrons move 
from 2s2p to 2s 2 at the very beginning of the decay of 
2s2p. The laser duration allows the Rabi oscillation to 
run a little more than half cycle, which bounces only a 
small fraction of the population back to 2s2p. The rela- 
tively few electrons returning to 2s2p have changed the 
phase by 7r through this Rabi flopping, thus reversing the 
sign of q in the resonance profile (See Ref. [l9j). When to 
increases, in Fig. UOf a). an additional ridge moves from 
the lower energy to the original peak position, stretching 
from 35.4 eV at to = 15 fs to 35.5 eV at t = 50 fs. This 
suggests that the laser divides the autoionization into two 
periods of time; the part prior to the laser is responsible 
for the regular Fano profile as shown by the "no laser" 
curve, and the part after the laser is responsible for the 
"inverse" shape that we have just discussed. Different 
delay times determine the fractions of these two parts, 
and the ridge forms at different energies as the interfer- 
ence pattern. Figure [TUT d) clearly shows the shifting of 
the ridge while the inverse peak stays at 35.55 eV. 

The 2s 2 profile in Fig. [TTT a) has a very gentle and 
monotonic attenuation along to, while the 2p 2 profile in 
Fig. Eta) features a "gap" in to where the profile drops to 
zero altogether before it regains some strength as to in- 
creases. This gap, as explained in Sec. IIII A 11 originates 
from the IR ionization; however in the A^ = 540 nm 
case, this mechanism is missing since the ionization rate 
for 2s 2 is low. Rabi oscillation and autoionization are 
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FIG. 11: (Color online) As Fig. [Tolbut for the 2s 2 resonance. 



the only influences on the evolution of 2s 2 . Since the 
Rabi oscillation runs only a little more than half cycle, it 
can be viewed as a "one way route" for electrons from 
2s2p to 2s 2 . Once arriving at 2s 2 , the electrons au- 
toionize quickly to form the 2s 2 profile. As to increases, 
more electrons will autoionize from 2s2p and less will be 
brought by the laser to 2s 2 , and the 2s 2 resonance be- 
comes weaker. 

The choice of 2s 2 as the second coupled state is opti- 
mal for practical measurement issues, such as the signal 
strength and energy resolution. One must remember that 
the 2s 2 resonance is on top of the lsep background if the 
total electron energies are measured. The maximum sig- 
nal intensity of the 2s 2 profile is nearly 10~ 8 eV^ 1 , and 
the nearby lsep background is 7.5 x 10~ 9 eV -1 , i.e., when 
the two are layered together, the resonance signal is high 
enough to stand out from the background. Furthermore, 
its resonance width Ti, = 0.125 cV is larger than that 
of 2s2p, making its visibility better even with the same 
energy resolution used in Rcf. 



24| . Considering both the 



height and the width in spectrum, detecting the 2s res- 
onance vs delay should be feasible. 

The complicatedly structured profiles seen in 
Figs. fTPTb-d) are in contrast to the Autler-Townes 
doublet in long-pulse cases. The latter features the 
splitting proportional to the field strength, and each 
peak in the doublet does not differ from the original 
single peak in shape. In order to clarify the difference 
between long and short pulses in a systematic way, in 
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FIG. 12: (Color online) Photoelectron profile of the 2s2p res- 
onance with a 540 nm, 1 ps laser pulse. The original Fano 
line-shape is plotted in the gray solid curve. 



Fig. Q21 we show what the 2s2p profile should look 
like if the laser is 1 ps long, with all other parameters 
unchanged. A 1 ps laser pulse is much longer than our 
atomic time scale and equivalent to a stable AC field, 
where time delay cannot be defined. The separation 
of the splitting is about the Rabi frequency 0.27 eV. 
Nevertheless, in a standard EIT setup, the \b) state is a 
bound state, which is not the case in the present system. 
As a result, in our example, the split peaks are fainter 
because of the finite lifetime of 2s 2 . 

If the laser duration shrinks from 1 ps to 50 fs, the mag- 
nitude of the splitting will be a function of the time delay. 
However, since the laser is still significantly longer than 
the 17 fs decay lifetime in concern, the delay-dependence 
is more relevant to the overlap between the pulses rather 
than the autoionization. This scenario has been studied 
with transient absorption spectroscopy in helium where 
the XUV is 30 fs and the IR is 42 fs The study 

showed that the splitting was tuned to maximum by over- 
lapping the pulses, and it disappeared when the pulses 
were totally separate. The system can be viewed as a 
dressed atom where the dressing condition can be slowly 
turned on or off by changing the time delay. 

As the comparison indicates, the mechanisms for the 
short and long pulses are indeed very different. For a 
long dressing pulse, the field enters as the coupling (off- 
diagonal) terms in the Hamiltonian of the system. When 
the Hamiltonian is diagonalized, the energy levels are 
shifted, which are then represented by the doublet. How- 
ever, for a short dressing pulse, the coupling strength 
changes quickly; the Rabi frequency is not well-defined, 
and the mechanism for the long pulse breaks down. In 
other words, the dressing field influences only a short 
temporal segment out of the whole autoionization, while 
most of the time, the resonance profile evolves without 
laser and aims at the single-peak Fano profile rather than 
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the doublet. 



IV. SUMMARY AND CONCLUSIONS 

A model has been developed to study the autoioniza- 
tion dynamics in a laser-dressed helium. An attosecond 
XUV excites the 2s2p resonance in a time-delayed 780 nm 
IR pulse. The IR can couple the 2s2p and 2p 2 states and 
ionize the two states. The photoclcctron energy spec- 
tra for different time delays are calculated and compared 
with the experiment t 24f . While the experimental energy 
resolution was not good enough to observe the resonance 
shape in detail, good agreement for the resonance peak 
intensity vs time delay between our model and the ex- 
periment has been achieved. The decay lifetime of 2s2p 
can be retrieved by this result. Because of the strong IR 
ionization, the coupling with 2p 2 is to open an efficient 
pathway where the 2s2p resonance can be depleted, i.e., 
the IR field modifies the profile with an overall depres- 
sion without changing the spectral shape, which is totally 
different from the typical three-level systems. 

To reduce the effect of ionization by the IR laser, we 
change the laser wavelength to 540 nm and consider the 



coupling of the 2s2p and 2s 2 states. The 2s 2 state has 
a larger binding energy where its ionization is negligi- 
ble. In this case, a complicated pattern in the 2s2p res- 
onance shape vs time delay has been found. The result 
is interpreted by the Rabi oscillation between the two 
autoionizing states whose cycles are confined by the 9 fs 
laser. In order to make connection of this result to the 
traditional dressed-atoms, wc change the laser duration 
to 1 ps and recover the Autler-Townes splitting in the 
spectrum, which also clarifies the difference between the 
short- and long-dressing fields. In order words, in the 
presence of an intense IR, autoionization dynamics can 
be changed significantly. The possibilities of such ma- 
nipulations are tremendous, and what one can gain from 
such experiments remains to be further explored. 
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